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Abstract 



We motivate lattice QCD studies of the parity-violating pion-nucleon cou- 
pling constant and extend flavor-conserving hadronic parity-violation from 
QCD to partially-quenched QCD. The parity-violating pion-nucleon coupling 
and the anapole form factor (and moment) of the proton are computed to 
one-loop order in the partially-quenched chiral expansion. For the parity- 
violating pion-nucleon interaction, we include the contributions from total 
derivative operators necessary to match the kinematics that will be used in 
lattice simulations. 
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I. INTRODUCTION 



While flavor- changing parity-violating (PV) interactions are well understood theoreti- 
cally and a great deal of precise data exists, knowledge of flavor-conserving parity-violation 
is rather primitive. Flavor-conserving parity-violation continues to be an area of intense 
investigation in the nuclear physics community. Its study is presently serving both to un- 
cover the structure of the nucleon in electron-scattering experiments such as SAMPLE [Q, 
and to determine PV flavor-conserving couplings between pions and nucleons [0,0. The 
problems that are encountered in this sector are both experimental and theoretical. On the 
experimental side, the PV signals, unlike those in flavor-changing processes, appear as small 
deviations in either a strong or an electromagnetic process, such as PV in ep — * ep, or in 
the circular polarization of the 7-ray emitted in ^^F* — F7 The current situation is 
somewhat confused by the fact that measurements of parity-violation in atoms and nuclei do 
not give rise to a consistent set of couplings between hadrons However, it is important 
to keep in mind that many of the "experimental" determinations of these couplings require 
theoretical inputs with varying degrees of reliability. Recently, it has been reemphasized that 
measurements of PV observables in the single-nucleon sector would significantly ameliorate 
the situation by eliminating many-body uncertainties 1^,^. Despite the inherent difficulty 
of such experiments, there are ongoing efforts to measure PV processes in systems with only 
one or two nucleons, such as the angular-asymmetry in np d'y Such measurements 

should provide a reliable determination of the leading-order (LO), momentum- independent 
weak itNN coupling constant, h^^j^j^. 

On the theoretical side, despite heroic efforts to model [jl^^ini] hadronic matrix elements 



of the four-quark operators that appear in the low-energy effective theory of the standard 
model, there are no reliable calculations of the PV couplings between hadrons. A first 
principles calculation of /ii-TVAf lattice QCD would therefore be extremely welcome. This 
would require a lattice QCD simulation of a correlator with three hadronic sources interacting 
via a four-quark operator. Unfortunately, chiral symmetry does not allow one to relate the 
irNN correlator to a correlator without the pion. On the bright side, the structure of the 
four-quark weak Hamiltonian requires a fiavor change in the nucleon and therefore there are 
no disconnected diagrams to be computed on the lattice. 

A standard difficulty in reliably computing h^^lfj^ on the lattice is that present computa- 
tional limitations necessitate the use of lattice quark masses that are significantly larger than 
those of nature. In order to make a connection between lattice calculations in the foreseeable 
future and nature, an extrapolation in the quark masses is required. Chiral perturbation 
theory (xPT) provides a systematic description of low-energy QCD near the chiral limit, 
and this technique has been extended to describe both quenched QCD (QQCD) |T^^16| and 
partially-quenched QCD (PQQCD) The hope is that future lattice simulations can 



be performed with quark masses that are small enough to guarantee a convergent chiral 
expansion, thus allowing a meaningful extrapolation to the quark masses of nature. Re- 



cently, it has been shown how to include the low- lying octet and decuplet of baryons |T9,20 



into partially-quenched chiral perturbation theory (PQxPT) |]T^,|TB|. In this work we will 
present /i^-tvat and the anapole moment of the proton at one-loop level in SU {4:\2)l0SU (4|2)ij 
PQxPT with two non-degenerate light quarks. Given that lattice simulations have on-shell 
particles in both the initial and final states, extraction of /i^^fAf requires the injection of 
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energy at the weak vertex. This energy non- conservation modifies the effective field theory 
description of the three-point function by requiring the inclusion of operators that are total 
derivatives. We present a calculation of /i^-TVAr ^^e one-loop level both for the case where 
energy is conserved and the pion is off its mass-shell with = 0, and for lattice kinematics 
where energy is injected at the weak vertex. 



II. HADRONIC PARITY VIOLATING INTERACTIONS 



The four-quark operators that contribute to flavor- conserving, low-energy PV processes 
can be classified by how they act in isospin space, AJ = 0, 1, 2. Their QCD evolution from 



the weak-scale down to the strong interaction scale has been computed previously ||10| , |21| , |22 
In this work we focus on the A/ = 1 interaction, as it is only in this channel that the 
LO operator in the chiral expansion is momentum independent. In QCD, the effective 
Lagrange density for A/ 
quark operators 



1 interactions at the quark- level is P^J2B| ], including strange 



C 



A/=l 



Gf sin 6„ 



V2 3 ^ 

where 6^ is the weak mixing angle, and where the four-quark operators are 



(1) 



9i 



O2 



Tl^lbqp flf^r^qa 
-- s"Ysf3q^-ff,-f5T^qa 
-- s"Yl5Si3q^-ft,T^qa 



(2) 



and their coefficients at the chiral symmetry breaking scale are [22 



C«(A,) 
5«(A,) 



1.10 , 0.068 , 0.234 , -0.697 ) 
5.61 , -1.90 , 4.74 , -2.67 ) 



(3) 



where we have neglected inter-generational mixing through the Cabibbo angle. In the limit 
of vanishing 6^^ and degenerate quark weak-isospin doublets, the standard model possesses a 
global symmetry that forbids A J = 1,2 parity- violation This symmetry is refiected in 
the relative magnitudes of the S^^^ versus the C^^\ If strange quarks were to play no role in 
the structure of the nucleon, two-flavor lattice QCD simulations would provide an accurate 
calculation of the matrix elements of ^i,,,4, but to the extent to which strange matrix elements 
are non- zero, the two-flavor simulations provide only part of the PV matrix elements. The 
role that strange quarks may play in these matrix elements has been considered in various 
unjustifled approaches, e.g. Ref. |T^, p2| and [0], and the strange quark is found to make 



a sizable contribution p4[ to the PV momentum independent interaction, /i^^^^r- However, 
the computation of these matrix elements in two-flavor simulations will be a vital flrst-step 
towards a rigorous understanding of flavor conserving, hadronic parity-violation. 

In order to match onto a chiral Lagrange density that describes low-energy PV processes, 
it is convenient to introduce the objects X£ j^, deflned as 
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Xa 
R 



(4) 



where 



C = exp 



. f , 



M 



7v^ 



TT ' 



(5) 



and where ^ — > L^f/^ = UC,R^, and — > UXI^U'^ under chiral transformations. In our 
convention, the pion decay constant, / = 132 MeV. At LO in the chiral expansion, AJ = 1 
PV couphngs between the pions, the nucleons and the A's are described by the Lagrange 
density @ Q, 



wk 



(1) / 

ttNN T 



N 



7,(1) 



- 



^(1) / Triabc,)! 
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T 



xl - xl 



Tr 



abdjfj. 



,(1) 



73 



■ttAA^Om A + ^ttAA -X— M A 



+ h.c. 



(6) 



III. EXTENSION TO PARTIALLY-QUENCHED QCD 

The extension of the PV flavor- conserving weak operators from QCD to PQQCD is anal- 
ogous to the extension of strangeness-changing operators, as discussed in detail in Ref. [26|. 



It is convenient to rewrite the operators in eq. (0) in terms of operators with well-defined 
properties under chiral transformations. Neglecting the strange-quark operators one finds, 

h = Tr'Y{l - I5)qa fl.il + 75)g/3 , h = rr'ril + l,)qa q^l,{l + l,)qp , (7) 

and operators 65.. Og, corresponding to the other color-contraction. These objects transform 
as (3, 1), (1,3), (3, 1), (1,3), respectively under SU{2)l ® SU{2)r chiral transformations. 
The most general extension of this operator set to PQQCD is quite messy, as one ends 
up dealing with objects of the form (A, A) of SU{4:\2)l ® SU{4:\2)r, where A denotes the 
adjoint representation of SU{4\2). Therefore, for simplicity reasons alone, we will only 
consider extensions of the form (A, 1) and (1, A) In this case, we can define the weak 
flavor matrix 



-3 



T 



-1 , hj , hi , hj , hi ) , (8) 



where hj and hi are arbitrary weak charges. This matrix is supertraceless and its matrix 
elements reduce to those of QCD when the sea-quark (j and /) masses become degenerate 



^Ref. uses similar notation with h^^^ = /i^ and h^J^^ = h^. 
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The flavor structure can be further extended when lattice simulations require. 
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with the valence-quark masses. The quark mass matrix in PQQCD is given by mq = 
diag(m„, m^, m^, m;, m^i). Further, in analogy with QCD we define 

(PQ)Xl = e , ^^^^X^ = ^ . (9) 

The LO Lagrange density describing the low-energy PV interactions in PQ^PT is 



wk 



-Wij (bb 

f 



{PQ)Xl 



B 



(10) 



where B, which contains the nucleons, and 7^, which contains the A-resonances, transform 
in the 70 and 44 representations of S'f/(4|2), respectively pOf. The contraction of fiavor 
indices, ( ), can be found in Ref. [|T^. At tree-level one can make the identification 



h 



(1) 

rNN 



Wo 



h 



(1) 

rAA 



(11) 



The LO weak coupling between the nucleons and A-resonances involves one derivative and 
hence is higher order in the chiral expansion. 

The Lagrange density describing the parity-conserving interactions of the 70 and 44 
with the pseudo-Goldstone bosons at LO in the chiral expansion is [|T5|] , 



C = 2a (BS^'BA.) + 2(3 [BS^A^B 



+ 



2niT'S^A^% 



[ (t'^a^b) + {BA'Z 



(12) 



where is the covariant spin- vector |p7H29||, and A^ = ^ ^d^^"^ — ^^d^^ )• A compari- 



son with the LO interaction Lagrange density of QCD yields the tree-level identification ||20 

/3 



a 



4 1 

-^9A + 



1 



^91 - ^9A 



n 



9AA 



c 



-9an 



(13) 



where qa, 9aa and qan are the NN, A A and AA^ isovector axial couplings, respectively, 
and {9a+9i) is the isosinglet axial coupling. For details about this and other aspects of 
PQxPT, we refer the reader to Ref. pOl . 



IV. hllfj^ AT ONE LOOP WITH = 

In QCD, the LO pion-nucleon PV interaction is generated by the Lagrange density 
in eq. (^. At higher orders, the PV pion-nucleon interaction will receive contributions 
from one-loop diagrams and also from terms at next-to-LO (NLO) in the chiral expansion 
involving a single insertion of ruq, 

C = -h{ NiMT"", Xl - Xl}N - cJ-n(xI- Xl) N tr (uf) . (14) 
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Since the pion mass is much greater than the proton-neutron mass difference, one does not 
have all three particles on their mass-shells. Therefore, we compute the one-loop corrections 
to the weak vertex in both QCD and PQQCD. The one-loop calculation of the vertex has 
been performed previously in QCD p5[, and explicit computation in the isospin limit and 
with vanishing pion four- momentum, P^, gives 



+ 



"'ttNN 
1 



2m (ci + C2) 



h 



(1) 

TvNN 



20 



(15) 



where m = rriu = rrid in the isospin limit. The loop functions are defined as 
m^log (m^//i^), and = J(m^, A,yu) with 



J(m, A, /i) = (m^ - 2A2) log ( — j + 2AVA2 - log 



'A- VA2-m2 + 2e' 
A + VA2-m2 + ie^ 



(16) 



once 



and A is the A-nucleon mass splitting. Our results agree with those of Ref. |2^ 
differences in conventions are accounted for. 

In PQQCD the counterterms involving a single insertion of tuq that contribute to the 
PV pion-nucleon interaction at the one-loop level are 



f 



(^_'^{m+nj){nk+vr.) -g^i^ | {pQ)xI -(^Q) x|, , M+}IB. 



+ C2 



+ 



+ C6 B'^'^ ((^«)Xl -(^'3) X|)' str (A< 



+ C7 



^iPQ)xl _{PQ) Xl)lB,,r, str (^_ 



+ C8B'^'^i3.,,strff(^«)X?,-(^«)X^l 



(17) 



where M.+ = | (^^mg^^ + ^mg^^. The vertex at one-loop level can be written as 

1 



iPQ)r_ 

rt ' 



p + 



where the diagrams in Fig. ^ give 



V° + hj 7]^ + hiT]^ ) + c° + hj c' + hi d 



(18) 



p=-- {2Wi - W2 ) 

7° = ~g [ Lju + Ljd + Liu + Lid + Rr]u,riu + Rvd,ild ~ 2-Rr;„,r;d 



+ 



9i9a 



[ (Wi - 2W2) {2Lud - Ljd - L,u - Lid - L 
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6 



I 



I 
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FIG. 1. One-loop graphs that give contributions of the form ~ ttt-q log mg to the momentum 
independent parity-violating interaction h^^^^. A solid, thick-solid and dashed line denote a 70-nu- 
cleon, 44-resonance, and a meson, respectively. The solid-squares denote an axial coupling given in 
eq. i^T^), while the crossed circle denotes an insertion of the parity-violating pion-nucleon operators 
with coefficients 1^1,2,3 in cq. (IC). The solid circle denotes an insertion of the strong two-pion 
vertex from the nucleon kinetic energy term. 



+3IVi {Lju + Ljd + Liu + Lid — '^Lud) ] 



9l 



+ 9\n n ^3 ( Jdd + Juu + 4J„d + Jjd + J ju + Jld 



-I-2T 



2T 
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■nn,Vd 



c° = - 
d = 



{ruu + nid) - ( -2ci + 4c2 - C3 - C4 + 2c5 
= , 



Ji 



lu 



[rrij + mi)-{ -Ce + 2c7 



(19) 



where R^^y = H{L^, Ly, Lx) and %^y = H{Jx, Jy, Jx), with 



nab{A,B,C) 



B 
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(20) 



where the mass, mx, is given by rrix = \ + The expression in eq. ([T9|) col- 



lapses down to the isospin-symmetric QCD expression given in eq. (0) in the limit, 
rrij, mi, rriu, rrid — * m. Wavefunction renormalization can be performed for those particles on 
their mass-shell using (expressions for Wp and Wn are given in Ref. and in the isospin 
limit is given in Ref. [[T^ ), 

"U^TT = 2 [ —Lju — Ljd — Liu — Lid + 2-Rr;„,r;d "~ Rvn,riu ~ Lirn^rn ] 
'^P ~ 9 A {Lud + Luu + 2Lju + 2Liu + 'iR-qu^riu) 

+ gi9A (2L 

+ {Luu — 5Lud + 2L;m + 3Lirf + 2Lju + 3Ljd + B-Rry^.rju + QRvu,Vd + 3-Rj7d,r?d) 



+ gi9A {'^Ldd — Lud + Ljd + Ljrf + "iRri^^rii + SRri^^jj^^ 



9l 



{Ldd — 5Lud + '^Ljd + ^Lju + 2Lid + 3Liu + 'iR-qu,vu + ^Rvu,Vd + '^Rvd.nd) 

+ -^9ANi^Lud+ Jdd + '2'Jju + '2'Jlu + ■Jjd + Jld + '2%,^^riu+'^'^Vd,Vd~ ^'^Vu,Vd) 5 (21) 



by adding a contribution 



'^^"^^r, = p ^, . (22) 



V. h^^lj^ AT ONE LOOP WITH LATTICE KINEMATICS 

The analysis of the previous section will not be as useful in determining /i^^^^Ar from the 
lattice as one would naively assume Q While the = limit is natural to use from the 
viewpoint of a momentum and rrig expansion, the fact that lattice simulations only measure 
on-shell to on-shell amplitudes means that the pions and nucleons are on their mass-shells in 
both the initial and final states. The extraction of /i^^rAr from Nit requires an injection 
of energy at the PV weak vertex which can occur because the weak operator is inserted on 
one time-slice only. Therefore, we must include contributions from operators that are total 
derivatives, which usually vanish. Recently, chiral perturbation theory has been used to 
describe K tttt with the kinematics appropriate for a lattice determination of the matrix 
elements of the relevant four-quark operators, m'^** = m^** and m'^** = 2m}^^^, including 



^We are indebted to Steve Sharpe for making this point clear to us. 
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the necessary total derivative terms [^. In QCD, the LO Lagrange density describing PV 
interactions is given in eq. (^), while the Lagrange density at NLO is 



{NLO) 
wk 



h^}) -ivDN 



xl - xl 



N 



(23) 



where is the nucleon four- velocity. This is the leading contribution from a heavy baryon re- 
duction of iD'^N'jfj^ [ X| — Xf^ ] N. Given baryon number conservation, the total derivative 
gives a non-zero contribution from the energy and momentum injected by the X| — X|. inser- 



tion. Working in the frame where the initial state nucleon (proton) is at rest, 
the amplitude at NLO resulting from eq. (^ and eq. (pSf) is 



A- 



{mT\i J t 



"-wNN ^ "' 



E 

J 



:i,0,0,0). 



(24) 



where E is the energy injected by the weak vertex. In order to produce an on-shell nir final 
state, the injected energy must exceed E > m.,^ + Mn — Mp. Near threshold, where the 
final state neutron and pion are at rest and E = m.,^ + Mn — Mp, the contribution from 

my^, and is formally dominant over the loop 



the total-derivative operator, scales as 



corrections and counterterms of the previous section. 



I 71 ™ 
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1^ 
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FIG. 2. Tree-level contributions to the parity-violating vertex with lattice kinematics. The 
crossed circle denotes an insertion of the parity-violating pion-nucleon operators with coefficients 
Wi 2 in eq. (f?^. The solid square denotes an insertion of the energy-momentum dependent opera- 



tors with coefficients Wdi,D2 in eq. (21 ). The crossed-box denotes an insertion of energy-momentum 
at the weak vertex. 



In PQQCD, the lagrange density describing PV interactions at NLO is 



— I : — V ■ D [B 



4 



{PQ)X^ _(PQ) ^3 



(25) 



The one-loop calculation of the previous section, with minor modifications, provides the 
leading non-analytic contributions arising at N^LO in the chiral expansion. All other con- 
tributions are formally suppressed in the chiral limit. The matrix element at one- loop can 
be written as, keeping only leading non-analytic terms. 



P + 



IQn^f^ 



rf - -p[wp + Wn + w^] + hj f]^ + hi rf 



Up, (26) 
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(I ( ) 



(e) I 



(f) 



M,X M ' 5C 

N ^ N N / \ N 



(g) (h) 

FIG. 3. One-loop graphs that give contributions of the form ~ mQlogmq to the momentum 
independent parity-violating interaction h^)^j^. A solid, thick-solid and dashed line denote a 70-nu- 
cleon, 44-resonance, and a meson, respectively. The small solid squares denote an axial coupling 



given in eq.(li), while the crossed circle denotes an insertion of the parity-violating pion-nucleon 
operators with coefficients Wi,2,3 in eq. j^ldj). The solid circle denotes an insertion of the strong 
two-pion vertex from the nucleon kinetic energy term. The crossed-box denotes an insertion of 
energy-momentum at the weak vertex. Diagrams (a) to (e) contribute to vertex renormalization 
while diagrams (f) to (h) contribute to wavefunction renormalization. 



where the wavefunction renormahzation contributions, Wi, can be found in eq. (ZT). The 
diagrams in Fig. |] give, retaining only the leading non-analytic contributions and working 
in the frame where the initial state proton is at rest. 
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'ud 
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(-) 
Id 
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(27) 



where p, rf., rj^ and rj^ are given in eq. ([T9|). The functions Z^j^"* = L{mij, ±E,jj,) are 
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L{m,+E,fi) = -AE 



E\og 



m 



— V E^ — m? loe 



'-E - y/W^m^^^e 
— m? + ie ^ 



(28) 



Note that these functions are enhanced by a chiral logarithm compared with contributions 
from local counterterms in the chiral limit. When E = ±171^^, corresponding to the pro- 
duction of a nucleon and pion at rest, L{m,±m,^) = — 4m^ log (m^//i^). The additional 
non-analytic contributions in eq. (^) result from modifications to diagrams (c) and (d) in 
Fig. ^, with the other diagrams unchanged. 



In the limit where m 



'j,l,u,d 



m. 



this matrix element reduces down to that of QCD, 
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Un 



-h. 
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■kNN 



- h 



E 
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167r2/^ 



) 

NN 
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We have not given the numerous counterterms that are expected to appear at N^LO 
whose /i-dependence will exactly cancel that of the one-loop expressions in eq. (pGf). The 
counterterms are expected to make contributions that are smaller than those of the one-loop 
graphs when a renormalization scale of ~ is chosen, where is the chiral symmetry 
breaking scale. 

For an energy injection of > m,r + Mn — Mp, the amplitude develops an imaginary 
part due to the rescattering of the pion in the intermediate state (final-state interactions). 
Therefore, in accordance with the Maiani- Testa theorem [^, the simulations must be done 
at threshold, with E = m-„ + M„ — Mp. 



VI. THE ANAPOLE MOMENT AND FORM FACTOR OF THE NUCLEON 

An object that plays an important role in PV eH scattering is the anapole moment of 
H, where H denotes a generic hadron. The dominant PV contributions to eH scattering 
arise not only from tree-level Z^-exchange between H and the electron, but also from the 
exchange of a photon along with hadronic PV from the four-quark operators in eq. (|^). For 
the proton, there is a contribution to the matrix element of the electromagnetic current in 
the presence of hadronic PV interactions, of the form 

{p\jL{q)\p) = Mq') Up {S^q' -S-qq^)Up , (30) 
p 

where Ap{q'^) denotes the anapole form factor (the anapole moment is defined to be Ap{0)) of 
the proton. While such electromagnetic contributions vanish for on-shell photons, they give 
rise to local operators involving the proton and e (or whatever charged probe is involved 
in the process). In QCD, the anapole moment and form factor of the proton have been 
computed in xPT p3| , |32| -p5[] , and are dominated by the /i|rJVAr coupling (assuming that it is 



not anomalously small compared to estimates based on naive dimensional analysis). Further, 
the anapole moments of a few nuclei have been studied theoretically (e.g. Ref. @]). The 
proton anapole form factor at LO in xPT is found to be []33| ,^ 
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487r/ 



12 



20mf 



+ ... 



mi+ — — I tan ^ f — 



(31) 



where the — > hmit reproduces the anapole moment calculation of Ref . |32| . 

While lattice computations of the proton anapole form factor are of somewhat secondary 
importance compared to the computation of /il-TVAr, it would be extremely interesting to have 
a lattice QCD determination. To this end, we present expressions for the anapole moment 
of the proton in two-flavor PQQCD. The LO contribution to the proton anapole form factor, 
arising from the one-loop diagrams shown in Fig. is 
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FIG. 4. One-loop graphs that give the leading contribution to the anapole moment and form 
factor of the nucleon. A solid and dashed line denote a 70-nucleon and a meson, respectively. The 
solid-squares denote an axial coupling given in eq. ^1^), while the crossed circle denotes an insertion 
of the parity-violating pion-nucleon operators with coefficients Wi^2 'in- eg. 47^ 
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Wi 


1 2^jiu 


-A 



Hd 



12 



+ ^ (W, {2Aud - Add - 2Ai^ + Aid) + W2 {2A^d - ^hd - "^Am + 4^4;^)) 
+ ^ (W^i (2ij„ + - 2i„„ - + 2W2 + 2i,d - Ku - 2ind)) 
+ 1^ {W. (2i. + - 2i„. - i..) + 2H^. (i. + 2i. - i.. - 2i..)) 

] , (32) 

where we have used the short-hand A^ = Axiq"^) and we have used the electric-charge matrix 
in PQQCD, 

= diag +^ , -1 , , qi, qj, qi^ , (33) 

where qj and qi are arbitrary charges. The form factor in eq. (^) clearly reduces to the 
QCD expression in eq. (|3l|) in the appropriate limit. 

VII. CONCLUSIONS 

With a new generation of experiments designed to determine the PV nucleon-nucleon in- 
teraction, it is time for a lattice QCD calculation of the flavor-conserving PV coupling /i^^^^^. 
Near future lattice simulations will be partially-quenched with quark masses significantly 
larger than those of nature. Therefore we have presented the one-loop level expressions re- 
quired for the extraction of /i^-TVAf) both with vanishing pion four- momentum and with lattice 
kinematics. Furthermore, we have -with extreme optimism- given the one-loop expressions 
necessary to determine the proton anapole moment and form factor. 
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